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The effects of imperfect gate operations in implementation of Shor's prime factorization algorithm 
are investigated. The gate imperfections may be classified into three categories: the systematic error, 
the random error, and the one with combined errors. It is found that Shor's algorithm is robust 
against the systematic errors but is vulnerable to the random errors. Error threshold is given to the 
algorithm for a given number N to be factorized. 
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I. INTRODUCTION 

Shor's factorization algorithm JjJ is a very important quantum algorithm, through which one has demonstrated the 
power of quantum computers. It has greatly promoted the worldwide research in quantum computing over the past 
few years. In practice, however, quantum systems are subject to influence of environment, and in addition, quantum 
gate operations are often imperfect Environment influence on the system can cause decoherence of quantum 

states, and gate imperfection leads to errors in quantum computing. Thanks to Shor's another important work, in 
which he showed that quantum error correlation can be corrected g. With quantum error correction scheme, errors 
arising from both decoherence and imperfection can be corrected. 

There have been several works on the effects of decoherence on Shor's algorithm. Sun et al. discussed the effect 
of decoherence on the algorithm by modeling the environment ||. Palma studied the effects of both decoherence 
and gate imperfection in ion trap quantum computers Q . There have also been many other studies on the quantum 
algorithm 

The error correction scheme uses available resources. Thus it is important to study the robustness of the algorithm 
O h> itself so that one can strike a balance between the amount of quantum error correction and the amount of qubits 
available. In this paper, we investigate the effects of gate imperfection on the efficiency of Shor's factorization 
• algorithm. The results may guide us in practice to suppress deliberately those errors that influence the algorithm 

most sensitively. For those errors that do not affect the algorithm very much, we may ignore them as a good 
approximation. In addition, study of the robustness of algorithm to errors is important where one can not apply the 
. - - i quantum error correction at all, for instance, in cases that there are not enough qubits available. 

The paper is organized as follows. Section II is devoted to an outline of Shor's algorithm and different error's 
modes. In Section III, we present the results. Finally, a short summary is given in Section IV. 
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II. SHOR'S ALGORITHM AND ERROR'S MODES 



Shor's algorithm consists of the following steps: 

1) preparing a superposition of evenly distributed states 
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where q — 2 L and N 2 < q < 2N 2 with TV being the number to be factorized; 
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2) implementing y°modiV and putting the results into the 2nd register 

l^i> = ^Ei a >iy amodiV >; 

3) making a measument on the 2nd register; The state of the register is then 
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where j < 
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4) performing discrete Fourier transformation (DFT) on the first register |</> 3 ) = (j2 c f (c) |c}^ |z), where 



/ (c) = — ^ exp 

q 3=0 



2m(jr + 1) 



e i > exp 
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This term is nonzero only when c = k^, with fc = 0, 1, 2...r — 1, which correspond to the peaks of the distribution in 
the measured results, and thus this term becomes /(c) = -^e " . The Fourier transformation is important because 
it makes the state in the first register the same for all possible values in the 2nd register. The DFT is constructed by 

two basic gate operations: the single bit gate operation Aj = ^ J ~1 ) ' w ^ c ^ * s a ^ so caucd tnc Walsh-Hadmard 
transformation, and the 2-bits controlled rotation 
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with 9jk = 2^rj- The gate sequence for implementing DFT is 

{Aq-l)(B q _2q-lA q -2) ■ ■ ■ (B 0q _iB 0q _ 2 ■ ■ ■ B lA ). 

Errors can occur in both Aj and Bjk- Aj is actually a rotation about y-axis through | 
If the gate operation is not perfect, the rotation is not exactly §. In this case, Aj is a rotation of J + 25 



M*) = 



cos((5) — sin((5) — (sin(<5) + cos(<5)) 
V Sm (^) + cos((5) cos(5) — sin(<5) 



If 5 is very small, we have: 
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Similarly, errors in Bjk can be written as 

/l 00 \ 
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Bjk ~ 1 
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With these errors, the DFT becomes 
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(1 + S' c )\c) = ± J2 e^ +5 ^(l + S' c )\~c), 



(i) 
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where 5 C and S' c denote the error of Aj and Bjk, respectively. 

Let us assume the following error modes: 1) systematic errors, where S c or S' c in (1) can only have systematic errors 
(EMi); 2) random errors (EM2), for which we assume that <5 C or S' c can only be random errors of the Gaussian or the 
uniform type; 3) coexistence of both systematic and random errors (EM 3 ). In the next section, we shall present the 
results of numerical simulations and discuss the effects of imperfect gate operation on the DFT algorithm, and thus 
on the Shor's algorithm. 



III. INFLUENCE OF IMPERFECT GATE OPERATIONS 



We first discuss the influence of imperfect gate operations in the initial preparation 

4-iA,_ 2 ...4)|0...0) - ^(|0) + |1) + MO) - |1») ® (|0) + |1) + fcflO) - |1))) ® . . . ® (|0) + |1> + S n {\0) |1») 

= -^rEi li 2...»n=ol ili2 - i "> + ^rELi i nE! 1 , 2 ..,„,o(l'i-»K-iOis+i..i»> - \h..i R - 1 li R+1 ..i n ) 

If the errors are systematic, for instance, caused by the inaccurate calibration of the rotations, then Si = 82 = ■ ■ ■ = 
S n = 5. In this case, we can write the 2nd term as 

I 1 
= -7=6 ^ (2s-n)|iii 2 -in>, 
v2 i 1 i 2 ...i„=a 

where s stands for the number of l's, and 2s — n = s — (n — s) is the difference in the number of l's and O's. Thus 
the results after the first procedure is 
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This implies that after the procedure, the amplitude of each state is no longer equal, but have slight difference. 
Combining the effect in the initialization and in the DFT, we have 

(1 + «y )(i + s c y^ +5 '^ = (i + <5'V (2 ^ +5> , 

where <5" = 5 C + 5 a . In the DFT, we have 

q c=0 j=0 

where we have rewrite 5" as Sj here. Let P c denote the probability of getting the state |c) after we perform a 
measurement, we have 



p < = 4EE( 1 + ^)( 1 + ^) e4(2f£+C)(mr+ ° xe " ( ^ +dl)( ' £r+ ' ) 

^ m=0 fe=0 

= 4 E + + **) cos[— r(m - fc) + (mr + OC - + 0«fc] (3) 

^ m k ^ 

From Eq. (3), we find that after the last measurement, each state can be extracted with a probability which is 
nonzero, and the offset I can't be eliminated. 

Eq. (3) is very complicated, so we will make some predigestions to discuss different error modes for convenience. 
Generally speaking, the influence of exponential error Sj is more remarkable than Sj, so we can omit the error Sj, thus 

DFT 9 |^> = f ElZli:Le i{ ^ +5 ^ jr+l) \c) • 
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A. Case 1 



If only systematic errors (EMi) are considered, namely all the Sj's are equal, then /(c) can be given analytically 
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The relative probability of finding c is 

and if c = k-, then 
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It can be easily seen that lim^o Pc = h, which is just the case that no error is considered. 

When S takes certain values, say, S — ^(k — ^)ir where k is an integer, then the summation in Eq. (|j) is on longer 

valid. In our simulation, 5 does not take these values. Here we consider the case where q = 2 7 = 128 and r = 4. For 
comparisons, we have drawn the relative probability for obtaining state c in Fig.l. for this given example. We have 
found the following results: 

(i) When 5 is small, the errors do hardly influence the final result, for instance when c = k-, then 

, rsin 2 (%) 1 
hm P c = hm t. — | — = - . 

<5_o <5^o q 2 sur(4p) r 

The probability distribution is almost identical to those without errors. 

(ii) Let us increase (5 gradually, from Fig. 2, we see that a gradual change in the probability distribution takes place. 
(Here, we again consider the relative probabilities) When S is increased to certain values, the positions of peaks change 
greatly. For instance at S = 0.05, there appears a peak at c=127, whereas it is P c = when no systematic errors are 
present. In general, the influence of systematic errors on the algorithm is a shift of the peak positions. This influences 
the final results directly. 

B. Case 2 

When both random errors and systematic errors are present, we add random errors to the simulation. To see the 
effect of different mode of random errors, we use two random number generators. One is the Gaussian mode and the 
other is the uniform mode. In this case, the error has the form 6 = Sq + s, where i5o is the systematic error, s has a 
probability distribution with respect to c, depending on the uniform or the Gaussian distribution. When So = 0, we 
have only random errors which is our error mode 2. When So ^ 0, we have error mode 3. For the uniform distribution, 
s ~ ±s ma x x u(0, 1) where u(0, 1) is evenly distributed in [0,1]. s max indicates the maximum deviation from Sq. For 
Gaussian distribution, s ~ N(0,ao). Through the figure, we see the following: 

(1) When only random errors are present (Sq — 0), the peak positions are not affected by these random errors. 
However, different random error modes cause similar results. The results for uniform random error mode are shown 
in Fig. 3. For the uniform distribution error mode, with increasing S max , the final probability distribution of the 
final results become irregular. In particular, when S ma x is very large, all the patterns are destroyed and is hardly 
recognizable. Many unexpected small peaks appear. For the Gaussian distribution error mode, as shown in Fig. 4, 
the influence of the error is more serious. This is because in Gaussian distribution, there is no cut-off of errors. Large 
errors can occur although their probability is small. The influence of txo on the final results is also sensitive, because 
it determines the shape of the distribution. When <7o increases, the final probability distribution becomes very messy. 
A small change in can cause a big change in the final results. 

(2) When 8$ ^ 0, which corresponds to error mode 3, the effect is seen as to shift the positions of the peaks in addition 
to the influences of the random errors. 
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IV. SUMMARY 



To summarize, we have analyzed the errors in Shor's factorization algorithm. It has been seen that the effect of the 
systematic errors is to shift the positions of the peaks, whereas the random errors change the shape of the probability 
distribution. For systematic errors, the shape of the distribution of the final results is hardly destroyed, though 
displaced. We can still use the result with several trial guesses to obtain the right results because the peak positions 
are shifted only slightly. However, the random errors are detrimental to the algorithm and should be reduced as much 
as possible. It is different from the case with Grover's algorithm where systematic errors are disastrous while random 
errors are less harmful 191. 



[1] P.W. Shor, Proceedings of the 35th Annual Symposium on the Foundations of Computer Science, edited by S. Goldwasser 

(IEEE Computer Society Press, Los Alamitos, CA, 1994) p. 124. 
[2] A. Ekert and R. Jozsa, Rev. Mod. Phys. 68 (1996) 733. 
[3] W.G. Unruh, Phys. Rev A51 (1995) 992. 



[4] I. Chuang and R. laflamme, "Quantum error correction by codding" (1995) :niant-ph/9511003. 
[5] CP. Sun, H. Zhan and X.F. Liu, Phys. Rew. A58 (1998) 1810. 

[6] GM. Palma, K.A. Suominen and A.K. Ekert, Proc. R. Soc. London, A 452 (1996) 567. 

[7] R.P. Feynman, Int. J. Theo. Phys., 21 (1982) 467. 

[8] D. Deutsch, Proc. R. Soc. Land. A 400 (1985) 97. 

[9] L.K. Grover, Phys. Rev, Lett, 79 (1997) 325. 
[10] GL. Long, Y.S. Li, W.L. Zhang, C.C. Tu, Phys. Rev. A 61 (2000) 042305. 
[11] L.K. Grover, Phys. Rev. Lett, 80 (1998) 4329. 

Figure Captions: 

Fig.l. Relative probability for finding state c in the absence of errors. 

Fig. 2. The same as Fig.l. with systematic errors. In sub-figures (1), (2), (3), (4), S are 0.02, 0.03, 0.05 respectively. 
In sub-figure (4), the curve with solid circles (with higher peaks) is the result with 5 = 0.1, and the one without solid 
circles(with lower peaks) denotes the result with S = 0.33. 

Fig. 3. The same as Fig.l. with uniform random errors. In sub-figures (1), (2), (3), (4), s max are set to 0.01, 0.03, 
0.05, 0.1 respectively. 

Fig. 4. The same as Fig.l. with Gaussian random errors and systematic errors. In sub-figures (1), (2), and (3) r are 
set to 0.01, 0.03 and 0.05 respectively, and So = 0(without systematic errors). In sub-figure (4), both systematic and 
random Gaussian errors exist, where Sq — 0.33, t = 0.02. 
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